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Abstract 

Keeping nickel, cobalt and iron in mind, we investigate the origin of the itin- 
erant ferromagnetism. Recent experiments show that the systems should be 
in an intermediate coupling regime, where the band width and the interaction 
energy are of the same order of magnitudes. To treat such a situation, we gen- 
eralize the Gutzwiller approximation. In that, we take account of the effect of 
the band degeneracy and the Hund's-rule coupling in addition to the on-site 
repulsion. In generalizing the Gutzwiller approximation to the bands with 
degeneracy, we introduce the intuitive way to give the required expressions. 
After the discussion on nickel, the condition for the incomplete ferromag- 
netism, observed in cobalt and iron, is argued. Phase diagrams, which show 
the interplay between the band shape peculiarity and the Hund's-rule cou- 
pling, are given. It is found that for the 3d-transition metal systems, both of 
the Hund's-rule coupling and the special feature of the density of states are 
necessary to explain the itinerant ferromagnetism. 



* E-mail : okabe@ton . scphys . kyoto-u . ac . j p 



1 



I. INTRODUCTION 



Recently, studies on the electron correlation are extensive. While this is undoubtedly 
due to the high Tc materials in addition to the heavy fermion compounds and the related 
topics, the first motivation to elucidate the nature of the highly correlated system was to 
clarify the origin of the itinerant ferromagnetism. Our purpose of this paper is to study 
on the origin of the itinerant ferromagnetism in nickel, cobalt and iron. 

On the one side, there exist a lot of works which treat the purely mathematical mod- 
els 1^, since it is clear now that the effect of the electron correlation plays an important 
role. Although these works contain the interesting insights into the mechanism for the fer- 
romagnetism, it may give the unrealistic impression. On the other side, we see the fact 
that the band calculation using the local density approximation, which includes roughly 
the electron correlation effects, explains well the observed value of the spontaneous mag- 
netization etc. for the metallic ferromagnets nickel, cobalt and iron. Still, we should note 
that we cannot identify what mechanism plays a central role for ferromagnetism in terms of 
the band calculation. Therefore, to elucidate the physics of the ferromagnetism, we should 
start with the model which still describes closely the real materials. We take into account 
of the effect of the band degeneracy and the Hund's-rule coupling. As the interaction is 
regarded as the on-site type in our model, the electron correlation included in our calcula- 
tion is thought to be included partially also in the band calculation. Therefore, our work 
can explain the successful result of the ground state property derived with use of the local 
density approximation. 

Although there exist the prominent model calculations since several decades ago, most 
of them assumed the strong coupling interaction. For example, Kanamori, in his famous 
work argued on the magnetism of nickel, where he considered the effect of the band 
degeneracy. His conclusion was that we can neglect the band degeneracy and we are allowed 
to treat a single band model. This was due to his assumption J <^ /, where J and / are 
respectively the exchange integral and the interaction energy between holes with the an- 
tiparallel spins. He also assumed / ~ 2W, where W is the band width. Recent experimental 
findings say that these are not true but that the Hund's-rule coupling is effective and 
that the strength of the interaction is not large compared with the band width. 

Thus we should retreat the problem of the itinerant ferromagnetism which is in the 
intermediate coupling regime of the correlated carriers in the degenerate bands. This is one 
of the motivation of our work. 

Since the condition derived by Kanamori is based on the modified Stoner condition. 
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p{eF)Ieff > 1, his conclusion does not necessarily show that nickel is completely ferromag- 
netic in the ground state. (It is not enough to prove the instability of the paramagnetic 
state.) Therefore, we should reinvestigate this point: the ground state energy over the 
whole range of the magnetization should be calculated to show the global stability of the 
complete ferromagnetic state. 

While Kanamori limited his argument to the case with dilute carriers, especially for 
nickel, we argue also on the incomplete ferromagnetism which is observed in such as iron 
and cobalt. Since they have more than one holes per site, we should adopt the other approx- 
imating scheme than the t-matrix theory by Kanamori. Extension of the approximation to 
discuss more generally the degenerate bands is the other purpose of our paper. 

Another aspect of our work is the reinvestigation of the old ideas of Van Vleck and 
Herring [^. 

Some decades ago. Van Vleck argued that the effect of the Hund's-rule coupling would 
be large enough to explain the itinerant ferromagnetism, even for the metallic nickel in which 
the hole number is so small as 0.6 per site. According to him, the polarity fluctuation plays 
an important role, that is, more than one carriers can be on a site at a time to gain the kinetic 
energy, and thereby the Hund's-rule coupling becomes relevant. Thus, for the Hund's-rule 
coupling to be effective, the inter-particle repulsion should be reduced by the screening 
effect enough to allow the nearly free carrier migration. Herring conjectured that this 
actually would be the case for nickel, and that the band degeneracy with the Hund's-rule 
coupling should be crucial for the ferromagnetism in nickel. (His conjecture with respect to 
the effectiveness of the screening was true in the light of the recent experiments.) In other 
words, for the "id transition metal system, the ferromagnetism would not be realized without 
the band degeneracy and the Hund's-rule coupling. We investigate on this point. 

Here in this paper, as an approximating scheme among others, we adopt the Gutzwiller 
approximation 0] which is suitable for our purpose: the procedure is physically sound and 
transparent in addition to the fact that it is unrestricted by the carrier density or the re- 
pulsion strength, although it is generally hard to improve on the approximation further. 
It is expected that the Gutzwiller approximation is adequate to estimate the ground state 
energy of the correlated three dimensional systems unless the correlation is strong enough 
to drive the system toward the insulating phase. Therefore, we expect that in the range 
of our interest, the approximation is legitimate. There exist the works due to Chao and 
Gutzwiller in which they generalized the Gutzwiller approximation to the doubly degen- 
erate bands. However, their investigation was restricted to the case in the strong correlation 
regime, where the events that more than two carriers are on a single site are ignored. Our 
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expressions obtained below are the generalization of their works. Here we do not assume 
any restriction with respect to the strength of the interaction as well as the degree of the 
degeneracy. 

In the subsequent sections (§2 and §3), we generalize the Gutzwiller approximation to 
treat the degenerate bands. Then, after we investigate on the ferromagnetism in nickel (§4), 
we argue on the incomplete ferromagnetism (§5). We conclude that both of the Hund's-rule 
coupling and the special feature of the density of states are necessary to explain the itinerant 
ferromagnetism of the 3c?-transition metals. 



II. GUTZWILLER APPROXIMATION 

As an approximating method to investigate the single band Hubbard model, 
Gutzwiller p[ made use of the variational principle. To evaluate the expectation values 
with the state represented by the so-called Gutzwiller wavefunction >= Hjl ~ (1 ~ 
g)nii^nii]\'^o >, where |\l/o > denotes the Slater determinant and g{< 1) is the Gutzwiller 
variational parameter, he introduced the further approximation, which now is known as the 
Gutzwiller approximation. In this paper, we primarily focus on this method. After the sev- 



eral illuminating studies on the Gutzwiller approximation, |11,10] it has become clear that 
the approximation can be reproduced by a simple counting process. Although it is simple 
enough that it can be directly applicable for many models, it still involve the great amount 
of task if the model concerned is general. Therefore, in this section, before working on the 
actual variational calculation, we give the general argument to apply the approximation to 
the general situations and introduce the intuitive way to derive the expressions which are 
required. Our method is simple so that we do not have to go through the troublesome 
counting procedures. We use neither the explicit expression for the Gutzwiller wavefunction 
nor the Gutzwiller variational parameter g, but our final expression for the ground state 
energy, to be minimized, is expressed exclusively by the probabilities per site such as d, of 
the double occupancy. 

As an example, we start with the Hubbard model. 

H = Hkin + Hint = t 4<7ajf7 + UY^ rii^rii^^. (1) 

<i,j>,cr i 

With the interaction Hint, the probability that a site i is doubly occupied, which we denote 
as di, is reduced from ni^jUi^^, the value for U = 0. Here njo- designates the probability that 
the a spin electron is on the i-th site. (Or the mean number of carriers with the a spin.) 
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This reduction causes the decrease of the expectation value of the hopping part of the 
energy, < \1/ 1 if^ j„ | \1/ >. In the Gutzwiller approximation, this effect of the band width 
narrowing is taken into account only through the factor q„ (< 1, independent of i and 
j), by which < \l/|a|^ajo-|\I' > is multiplied. This factor q„ should be given in terms of d 
(independent of i) and the total energy E{d) is to be minimized with respect to d. 

E{d)=Y,q,{d)e, + Ud, (2) 

a 

where 

ea = t J2 < ^olalaaj^l^o >, (3) 

<i,j> 

is the kinetic energy of the uncorrelated state |\E'o > for U = 0. 

We begin with rewriting for < aj^aj^ >, obtained originally by Gutzwiller 0|, as 



— n + d\/n„ — d + \Jn_fj — dy/d^ 
nJl - rirr) 



+ v^j'^Vd'^iv^j'^v^^''^ + (4) 

(5) 



where n = rio- + n^fj is the number of electrons per site, = n^j — d and v = 1 — n + d 
denote the probabilities that a site is singly occupied by the a spin electron and that a 
site is vacant respectively. Suffices i and j represent the sites between which the hopping 
process takes place. The superscripts (6) and (a) are attached to indicate the probability 
amplitudes (square root of the probability) before and after the hopping process respectively. 
For example, the terms (^sj^''^ v^^''^)^(v^^*^ for q^ and {^^^^'^Vdt^'^)i{Vdt^^^^'''^)j 

for g| in eq. (§) represent the hopping process shown in Fig.l a) and b) respectively. 

The denominator takes just the same value as the numerator for d = n„n^„] n„{l — n„) 
might as well be written as (\/l — n^^''^ ^/n^^^^)i{y/n^^''^ y/1 — nj~°''')j. Thus g^- = 1 for [/ = 
as it should be. 

Here the point should be stressed: 

While our derivation of eq. (^) ( and eq. (|^) below) is intuitive and is not justified by 
itself, we can readily show that the final expressions derived after the lengthy calculations 
following the usual procedure of the Gutzwiller approximation can be cast into the form 
which can be interpreted as we propose here. In effect, the way we write down them (which 
is by no means self-evident at the outset) is noticed through the course of the involved 
calculations. Therefore, the approximating scheme itself does not contain any novelty but the 
original idea of Gutzwiller. However, it is a great convenience to know such an interpretation. 
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Hereafter we use d as the variational parameter and the other variational parameter g 
[r] in Gutzwiller's original papers) is not considered explicitly. It is easily seen that, in the 
Gutzwiller approximation, g is given by 

2 ^ {l-n + d)d 
{ua- - d){n-a - dy 

and the condition g <1 is equivalent to d < ricrn^cr, which is physically obvious. It is clear 
from the above derivation that is independent not only of the sites i and j but also of the 
type of operators, whether creation or annihilation. That is to say, if < \E'|a|o.a]^|\I' > (or < 
^laiffCtjcrl^ >), for example, takes non-zero value for U = 0, in the Gutzwiller approximation, 
this expectation value also is multiplied by the same factor when Hint is taken into 
account. [|l^] 

Now, it is straightforward to generalize the above procedure. As a result, the final 
expression for the total energy is given as 

h P>2{ii,-,ip} 

where u^^^li, ■ ■ ■ ,lp) denotes the probability per site that p electrons occupy the p states, 
{h, - ■ ■ ,lp} , with the other states being empty. (We use the indices Ip to discern the orbital as 
well as the spin states.) C^^\li, ■■■ ,lp) is the intra-site repulsion energy for the configuration 
u^P^li, ■ ■ ■ , Ip). The second summation sums over the different combinations of {Zi, ■ ■ ■ , Ip}. 
(Note that the Pauh principle requires 1^ 7^ in for m 7^ n.) e{l) (< 0) is the band energy for 
C^p) = (in case with no correlation), where the origin of the energy is taken such that the 
completely filled band gives 

e{l) = 0. (for the filled band) (8) 

Note that e is defined as a sum over occupied states of the kinetic energy e of each particle. 
(See eqs.(2.3), (2.12) and (2.14).) The effect of the electron correlation is included in the 
factor q{l), by which the band width is reduced. As above, the factor q{l) can be written 
down as 

/ \ 2 

1 



where n{li) is the average electron number in the /i state. Note that n{li) 7^ i^{h), but 

n{h) = iy{h)+ E E z/(p)(/i,/2,---,g- (10) 

p{>2) {h,-,lp} 
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The parameters {z/^^^} are determined by minimizing the energy E({z/(P)}), eq. (13), under 
the condition of the probabihty conservation, 

E E z.(^)(/i,---,g = i, (11) 

and the number conservation, eq. ([T0|), for given n{l). 

As for e{l), it is convenient to give it as a function of the carrier density. In principle, 
this is easily achieved: 

e{l) = p tpi{t)de. (12) 

n{l) = p pi{e)de. (13) 

If the density of states p/(e) is given, e{l) can be regarded as an implicit function of n{l) 
after eliminating e^, the Fermi energy. Here, as mentioned above (eq. (^) the origin of the 
energy should be given such that 

/oo 
epi{e)de = 0. (14) 
-oo 

Below we investigate on the ferromagnetism by calculating the energy as a function of the 
magnetic moment \n^ ~ "^J, provided that the shape of the density of states and the total 
carrier density n = + are given. 

In passing, to conclude this section, we give the simple argument to generalize the 
Gutzwiller approximation to the antiferromagnetic case. Here, for simplicity, we assume 
a simple two-sublattice (AB) antiferromagnetic structure. In the Hartree-Fock approxima- 
tion, the staggered magnetization m is given as 

-(A) = ^ E (15) 



The gap A is given by the minimization of the energy 

U 



where 



Ehf{A) = e{A) + ^{n^ - m^), (16) 



is the kinetic energy for the antiferromagnetic state. In the Gutzwiller approximation, the 
interaction term of eq. (|l^) should be diminished by the reduction of = (i^-t^)(^^-^) = 



7 



naUi) to d. As a result of this reduction, the kinetic term should be modified by a factor q. 
In the same manner as above (eq. (^), the factor q (independent of the spin component) 
may be written down as 

(Vl - n + dy/ua - d + y/ub - d\fd){^nb - dy/1 - n + d + Vdy/ua - d) 

Q= 1 ) (18) 

^Jriail -na)nb{l - rib) 

where Ua = ^^^^ and Ub = ^^^^^ are the number per site of the electrons with the spin 
component parallel and antiparallel to the staggered magnetization, respectively, {q = 1 for 
d = UaUb as expected.) In consequence, the variational parameters are A and and are 
determined by minimizing the energy, 

E{A,d) = qe{A) + Ud (19) 

Note that the above expressions are the same as given by Kotliar and Ruckenstein and 
Metzner \Ti\, who solved the long-standing problem of generalizing the Gutzwiller approxi- 



mation in a satisfactory manner to antiferromagnetism. 

To conclude this section, we repeatedly stress that our intuitive derivation of the factor 
q is based on the laborious but simple counting procedure. Our interpretation was given for 
convenience' sake hindsight. 



III. FORMULAE 

Neglecting the orbital dependence of the interaction energy, we parameterize the on-site 
interaction only by the two parameters / and I — J, which represent the intra-site repulsion 
between the two carriers with the antiparallel and the parallel spin components, respectively. 
The Hund's-rule coupling is taken into account as the exchange integral J. For the bands 
with the D-io\d degeneracy, we assume that the kinetic energy of the uncorrelated band can 
be written as De{^), which means that all of the D bands have the same structures. As for 
the filling, we assume no orbital ordering [TH] and regard each of D orbitals as equivalent. In 



this case, the expression for the total energy can be written down in terms of the parameter 
Up^q which denotes the probability for each of the configurations which has the p up-spin and 
q down-spin electrons on a site. (Compare with P{p,q), below eq. (^6]).) 

„ P+q=2D 

E{{i^}) = E + dCp- nCqVp^q {pql + {pC, + qC,){I - J)} , (20) 

'^=Ta p,<?(p+<?>2) 

where 
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and 



fp+q=2D 

Ql = "i /I _ "i N X! ■ D-lCq-l^Vp,qy/Vp^q-l . (22) 



D V-'- D \p=0,(j=l 

(We used the notation nCm = mhn-mv ' which counts the combination of the orbitals with 



M\{N-M)\ ' 

the same energy.) 

The conservation laws (eqs. (|10D and ([TT|)) require 



and 



i^l,0 = ^ ~ X/ -D-lC'p-l • oCql'p^q, (23) 
P>'?{p+'?>2) 

p+i7=2_D 

t'O.l = ^ — X! ■ D-lC*g-l^'p,g, (24) 

P,'?(p+g>2) 



p+g=2D 

1^0,0 = l- n+ (P + Q- ^)dCp ■ oCqUp^q, (25) 

P>9(p+'7>2) 



where = ?t,| + n^. Therefore, ehminating ui^, z/o,i and z/o,o by eq. (PHj), ( PH ) and 
we should minimize eqs. (EUI) with respect to z/p ^ with p + q >2. After elimination of Up g, 



the energy E' can be regarded as a function of n|(or n^). We then seek for the minimum of 
E{ni^) to study magnetism. 

We remark here that the probability P{p, q) that a site is occupied by the p up-spin and 
q down-spin electrons is given by 

P{P, q) = nCp ■ oCqUp^q. (26) 

In the case where C{p, q) = for all p and q, we can give P{p, q) explicitly by 

P{v.q) = DCp-oCq{^y {^)\ (27) 

which is the expression that Van Vleck used to estimate the degree of the polarity 
fluctuation. Then, it is easily seen that 

P>9>1 

as expected, (eq. (0.)) 
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IV. NICKEL 



Although we studied the itinerant ferromagnetism in nickel previously [|T6|, our consid- 
eration there was under the restriction that only up to two carriers can sit on the same 
site at most as in the work of Chao (In our notation above, only 1^2,0, ^1,1 and were 
taken into account as variational parameters.) While this is not a bad assumption in view 
of the smallness of the carrier density, 0.6 per site, we can improve upon in this respect 
in terms of the results of the previous section. The other reason for the reinvestigation for 
nickel here is that there we overestimated the value of the Hund's-rule coupling J. Here 
in this paper, as parameters for nickel, we use 0] n ~ 0.6eV, W ~ 4eV, / ~2.4eV and 
J ~1.2eV, i.e. , i = I /W ~ 0.6 and j = J/W ~ 0.3. Anyway, the result remains the same: 
in the ferromagnetic state, about 10% of the probability that a site has two carriers (holes) 
{Ep,q,p+q=2PiP, q) = ^(2, 0) + P(l, 1) + P(0, 2) ~ 0.1 for / = J = 0), is not so modified by 
including the effect of / and J. That is to say, the on-site interaction energies to be adopted 
in our calculation has already been screened by conduction electrons in such a manner that 
/ is reduced to be comparable with the band width W but keeping the effectiveness of J 
nearly intact. 0] 

To discuss the ferromagnetism of nickel and palladium, Kanamori |[l[] took account of the 
fact that the holes in the d band occupy states of the Brillouin zone close to the point X, 
where the symmetry of the band-edge state is X5. The six X5 states (two degenerate states 
at each of the three equivalent X points) are made up of the three atomic functions of the 
t2g class ( xy, yz and zx symmetries). Consequently, we regard the system as consisting 
of the triply degenerate bands, i.e., D = 3 as Kanamori. Then we have 13 variational 
parameters, Up^q with 2 < p + q < 2D. As for the kinetic energy, it is known that the 
density of states shows a significant peak near the top of the band at the Fermi energy 
for a hypothetical nonmagnetic nickel and also for palladium. Kanamori argued that the 
quantitative difference of the density of states at the Fermi energy is one of the reasons by 
which palladium remains paramagnetic while the ferromagnetism is realized in nickel. Since 
it is clear that such a specific feature of the density of states is important, we take this into 
account by the schematic density of states p as a function of the carrier number per state, 
n, as 



p{n) 



Ci l>n>2no 

(29) 

2no > n > 



r + (l-r)(i^)2 



The Fermi energy e as a function of n can then be given as 
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■n 1 



dn + C2. 



(30) 



pin) 



The normalization constant Ci is determined by the requirement that 



(31) 



by the definition of the band width W . And the origin of the energy C2 is given such that 



In Fig. 2 and Fig.3, we show p as a function of n and e, where we set = 0.1, W = 1 and 
= 2. Below we fix uq = 0.1 so that the Fermi energy of the hypothetical paramagnetic 
nickel is just at the peak, and use the relative height of the peak, r~^, as a parameter which 
characterizes the density of states of nickel. 

In Fig. 4 we show the phase diagram for j = J/W = 0.3 obtained with use of the 
Gutzwiller approximation, together with the results of the Hartree-Fock approximation. 
Compared with the latter, we see that the paramagnetic phase is stabilized remarkably. 
This refiects the fact that the one-body approximation should not be taken seriously with 
regard to the appearance of the ferromagnetism. We also see that the phase boundaries for 
the Gutzwiller approximation are almost independent of i = I/W. This is because of the 
fact that the strength of the effective interaction (which cannot exceed the band width) does 
not depend strongly on the bare one in the strong coupling regime. In a way, this can be 
taken as representing the fact that the peak at the density of states, the factor rather 
than the Hund's-rule coupling J, plays decisive role to realize ferromagnetism in nickel. 
However, we should note that it is due to the Hund's-rule coupling effect which makes the 
reasonable range of r^^ ~ 2 to be the complete ferromagnetic region. In Figs. 5 and 6, we 
show the phase diagrams for j = 0.15 and j = 0, respectively. Without the Hund's-rule 
coupling, or J = 0, the argument essentially renders back to the case for the non-degenerate 
band, for which it is known that ferromagnetism is hard to be realized in a realistic range of 
parameters. Note that the condition derived by Kanamori corresponds to the dotted curve, 
which is the boundary between the paramagnetic and the incomplete ferromagnetic state, 
since Kanamori's argument is based on the modified Stoner condition, p(ei?)/e// > 1. In this 
respect, his conclusion for nickel should be re- investigated because actual nickel is completely 




(32) 



according to eq. ([T^). The kinetic energy e is now given as a function of n by 




(33) 
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ferromagnetic. We see that, to obtain the complete ferromagnetic state without the Hund's- 
rule couphng, should be larger than 5 even for z ~ 1. It is, therefore, legitimate to say 
that the Hund's-rule coupling plays an decisive role for the complete ferromagnetism. 

We should also note on the continuous transition from paramagnetism to complete fer- 
romagnetism. Generally speaking, the boundary between the complete and incomplete fer- 
romagnetism is independent of that between the incomplete ferromagnetism and the para- 
magnetism. Therefore the two boundaries can cross at some strength of the interaction. 
Dashed curve in Fig. 4 represents the first order transition line around which the both of the 
complete ferromagnetic and the paramagnetic state can be locally stable. Such a boundary 
arises in the region of the strong interaction, that is to say, likely to appear without the 
Hund's-rule coupling nor the peculiarity of the density of states. 

In Fig. 7, we show the magnetic moment — nil as a function of n = + n| for i = 0.6, 
i — j = 0.3, rio = 0.1 and = 2.4. We see that the complete ferromagnetism appears for 
the values of n for which the Fermi surface lies in the peak portion, while above or below 
some definite values of n ( n^l.Ol and n^O.27) the state becomes unstable to result in the 
precipitous decrease of the magnetic moment. 

As for the factor q, since the electron correlation is not strong, the band width does 
not change drastically. For example, for i = 0.6, i — j = 0.3, uq = 0.1 and = 2.4, the 
factor q is 0.98 for the complete ferromagnetic state while q = 0.94 for the hypothetical 
paramagnetic state. Accordingly, the probabilities, eq. (^), obtained variationally do not 
modified appreciably from those for the uncorrelated values [|1^] , eq. (^) . If we approximate 
the exchange splitting A by I '^di^ '' I"t=0' which is valid if the carrier density n is small, we 
can compare the result with the Hartree-Fock and the Gutzwiller approximation. The result 
is Ahf = 3.65AGutzwmer for the same parameters as above. 



V. INCOMPLETE FERROMAGNETISM 

In this section, we investigate the incomplete ferromagnetism, which is observed in cobalt 
and iron. To the author's knowledge, the quantitative estimation of the parameters such as 
the Hund's-rule coupling are not known experimentally. Therefore, noting the fact that the 
strength of the Hund's-rule coupling would not drastically change from nickel to iron, we 
proceed to see what should be the condition which makes the complete ferromagnetic state 
unstable, while keeping the paramagnetic state also unstable. Here we treat the hypothetical 
metal of the bands with the five-fold degeneracy, D = 5, with the carrier density n = 2.0, 
keeping in mind the metallic cobalt and iron. As variational parameters, we use Up g with 
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p + q < 5, while we regard i'p,q=^ for p + g > 6. This is legitimate because the charge 
fluctuation involving more than 5 carriers at a single site is quite small for n = 2 even in 
the uncorrelated case. It is easily imagined that the Hund's-rule coupling would be more 
effective than the case of nickel because of the large number density of the carrier n > 1. 

On the other hand, the density of states at the Fermi level takes the relatively high value, 
as in nickel. Another fact which is apparent in the band calculation for cobalt and iron is 



the double peak structure of the density of states. |]18[ The second peak which is away from 
the Fermi level can be effective in the presence of the electron correlation, while, in the 
one-body approximation, the condition for the ferromagnetism can be written down only 
with the value at the Fermi surface. To take these into account, we use the model density 
of states as in the previous section: 



p{n) 



Ci 1 > n > 4no 

s + {l-l){:^Y 4no > n > 2no ^ (34) 

2^0 > ^ > 



no 



where no=0.2. (See Fig. 8 and Fig. 9.) The second peak at n = 3rio (for s"^ > 1) affects 
the origin of the energy, C2. (eq. (|T^ .) Consequently the presence of the second peak 
works against the ferromagnetism, or spreads the incomplete ferromagnetic phase space into 
the complete ferromagnetism. The phase diagrams are given in Figs. 10, 11 and 12. We 
see the effects of the Hund's-rule coupling and that of the second peak (s^^ > 1). Note 
that the Hund's-rule coupling (j ~ 0.17) should not be large compared with the case for 
nickel. This is partly expected by the fact that the parameter j is defined by j = J/W, 
and the band width W (~ 5eV) is larger for cobalt and iron than for nickel {W ~ 4eV). 
(The wavefunctions for (i-electrons spread more than those for Ni, thus, the values of J and 
/ for Fe and Co are more reduced than those for Ni.) We also see that the phase boundary 
is sensitive to the Hund's-rule coupling j, as expected by the large carrier density. The 
incomplete ferromagnetic region is enlarged by the effect of the second peak as well as by 
the Hund's-rule coupling. As for the reduction of the band width, the factor g^O.9 is given 
in the region of our interest. This is in accordance with the result for nickel: the electron 
correlation is not strong. 

To conclude, for the incomplete ferromagnetism to be realized in the realistic range of 
~ 2, it is necessary that the Hund's-rule coupling should be j ~ 0.17 (see Fig.13) and the 
shape of the density of states need not have the same feature, favorable for ferromagnetism. 
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as for nickel. Especially the double peak structure can make the incomplete ferromagnetic 
state stable. With respect to this point, our conclusion notes that the ferromagnetism in 
cobalt and iron should not be taken to be due to the high density of states at the Fermi 
level, rather it is ascribed to the effect of the Hund's-rule coupling. With the help of the 
itinerant carrier motion, the ferromagnetic correlation can be wide-spread and result in the 
long-range ordering. In these cases with more than one carriers per site, the mechanism due 
to the polarity fluctuation assisted by the Hund's-rule coupling gives the clear and natural 
explanation of the itinerant ferromagnetism. 

VI. SUMMARY AND DISCUSSIONS 

We generalized the Gutzwiller approximation to treat the realistic model of the itinerant 
ferromagnetic systems, including the band degeneracy and the Hund's-rule coupling. In 
deriving the required expression for the reduction factor q of the band width, we found 
that it can be interpreted intuitively that q is proportional to the square of the sum of the 
multiple of the probability amplitudes (square root of the probabilities), (eq. (^.) 

We see that the Hund's-rule coupling J is necessary to explain the itinerant ferromag- 
netism in the 3-(i transition metal. The phase diagrams for n = 2 showed that the bound- 
ary between the paramagnetic and ferromagnetic states is sensitive to the strength of the 
Hund's-rule coupling. Specifically, j = J/W should be around 0.17 for n = 2. In general, it 
is difficult to realize the incomplete ferromagnetic state as the ground state, because it needs 
to destabilize both of the paramagnetic and the complete ferromagnetic states. In effect, the 
energy difference between the paramagnetic and ferromagnetic states is quite small, so that 
any effect such as the peculiarity of the band density of states plays an important role to 
determine the magnetism of the ground state. As we saw above, the double peak structure 
can affect the magnetism by the many-body effect, even though the density of states at the 
Fermi surface is fixed. 

Concerning the reduction factor of the band width, our result showed that the factor q 
does not change appreciably from unity. This is partly because of the smallness of the inter- 
action strength. Moreover, we should note that our approximation is not accurate enough 
to be compared with the other method such as the diagram technique and the t-matrix the- 
ory. |]19[ In fact, our method cannot describe correctly the mass enhancement effect around 
the Fermi surface, for it is the dynamical effect due to the many-body correlation. Therefore, 
it would not be proper here to argue on the discrepancy with respect to the reduction of 
the band width between the experiments and the result of the band calculation. As it 
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is the difficult task to treat tfie electron correlation in the band with degeneracy, it is the 
future problem to improve upon the method itself to give more quantitative discussion. 
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1 : 


Examples of the hopping process. 
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2 : 


p as a function of n. eq. (4.1)) 
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3 : 


p as a function of e. 
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4 : 


The phase diagram for _7=0.3. 
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5 : 


The phase diagram for j=0.15. 
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6 : 


The phase diagram for j—0. 
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7 : 


Magnetic moment as a function of n. 


Fig. 


8 : 


p as a function of n. (eq. (5.1)) 


Fig. 


9 : 


p as a function of e. 
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10 : 


The phase diagram for s~^=l and j ~ 0.17. 


Fig. 


11 : 


The phase diagram for s"^ — r~^/1.2 and j — 0.17 


Fig. 


12 : 


The phase diagram for s~^ = r~^/1.2 and j = 0.18 


Fig. 


13 : 


The phase diagram for — 2, — r~^/1.2 
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